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Abstract:  To  support  an  effort  to  estimate  characteristics  of  surface-connected  damage  (i.e.,  cracks)  in 
aerospace  and  propulsion  structures,  a  computer-controlled  scanning  capability  that  acquires  experimental  eddy 
current  data  was  developed  based  on  an  Agilent  4294A  Impedance  Analyzer.  This  instrument  is  being  used  to 
assess  performance  of  different  computational  electromagnetic  codes  in  terms  of  determining  how  well  they 
predict  the  complex  response  from  simulated  flaws  in  conducting  plates.  The  particular  benchmark  problem  of 
interest  in  this  paper  measures  the  real  and  imaginary  impedance  due  to  Electrical-Discharge  Machined  (EDM) 
notches  in  aluminum  plates  that  are  approximately  13  mm  thick.  These  EDM  notches  are  5mm  deep  and  have 
the  same  surface  length  and  width.  A  range  of  notches  have  been  generated  for  this  study,  which  includes 
vertically  oriented  notches  and  angled  notches  of  10,  20  and  30  degrees  relative  to  the  surface  normal. 

Keywords:  flaw  characterization,  complex  impedance,  eddy  current,  volume  integral  method. 


1.  Introduction 

In  the  last  few  decades,  several  eddy  current  simulation  tools  have  been  developed  both  commercially  and  as 
research  grade  software  packages.  These  are  based  on  several  different  numerical  methods  such  as  the  finite 
difference  time  domain  method,  finite  element,  boundary  element,  volume  integral,  and  mesh  less  methods. 
This  paper  focuses  on  a  previous  experiment  conducted  by  Burke  [1].  This  is  a  simple  problem  that  can  be  used 
to  determine  the  accuracy  and  efficiency  of  different  numerical  methods.  The  two  codes  used  in  this  paper  are 
VIC  3D  [2],  a  volume  integral  code,  and  EC  SIM  [3-4],  which  uses  the  boundary  element  method  for  solution. 
Future  experiments  will  be  provided  to  the  wider  community  for  assessment  of  simulation  codes  based  on 
different  numerical  methods. 


2.  Background  /  Benchmark  Experiment 

This  experiment  uses  a  large  circular  air  cored  probe,  scanned  over  the  length  of  a  large  rectangular  notch  in 
a  plate  of  aluminum.  The  overall  experimental  setup  can  be  seen  in  Figure  1,  with  the  dimensions  labeled  and 
given  in  Table  1.  The  proportions  of  the  plate,  the  probe,  and  the  test  parameters  have  been  chosen  such  that  the 
outer  radius  of  the  probe  and  the  length  of  the  notch  are  on  the  same  order.  For  each  scan  position  over  the 
plate,  the  change  in  impedance  (AZ)  is  measured.  Three  plates  were  used  for  this  experiment,  each  of  a  different 
alloy  chemistry  and  conductivity,  which  differs  from  the  alloy  used  in  the  original  Burke  study.  The  three  alloys 
used  were  6061  T6,  7075  T6  and  2024  T3.  Each  square  plate  has  notches  placed  in  the  center.  Figure  2  shows 
the  testing  apparatus  with  the  probe  and  plate  in  place.  This  experiment  was  conducted  according  to  the 
procedure  published  by  Harrison  et  al  [5].  Other  methods  that  use  voltage  measuring  instruments  have  been 
proposed,  but  these  all  rely  on  a  transfer  function  approach  between  the  measured  impedance  and  voltage  output 
[6].  The  concern  is  that  error  in  the  model  can  essentially  be  absorbed  by  the  transfer  function  and  therefore  is 
not  suitable  for  model  validation.  It  should  be  noted  however,  that  once  proper  benchmark  experiments  are 
conducted,  such  transfer  functions  have  great  value  in  making  field  instruments  compatible  with  simulations  and 
model-based  inverse  methods.  Table  1  lists  the  conductivity  values  used  for  this  experiment. 
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F igure  1 .  Experimental  setup 


Inner  radius 

6.15mm 

Outer  radius 

124  mm 

Length 

6.15  mm 

Number  of  turns 

3790 

Lift-off 

Test  specimen 

0.88  mm 

Conductivity 

3.06  S/m 

Plate  thickness 

Defect 

12.22  mm 

Length 

12.60  mm 

Depth 

5.00  mm 

Width 

Other 

0  28  mm 

Frequency 

900  Hz 

Tablet.  Experimental  parameters 


Figure  2.  Testing  Apparatus  with  probe  and  plate  in  position 
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3.  Experimental  Procedure 


Preliminary  measurements  are  needed  for  later  post-processing  of  the  data.  This  includes  the  DC  values  of 
resistance  and  inductance,  as  well  as  multiple  scans  to  determine  if  repeatability  is  affected  by  temperature 
variations  or  other  experimental  factors.  The  DC  value  for  probe  resistance  was  found  using  a  4  point  analysis. 
The  DC  value  for  inductance  was  found  by  running  the  probe  over  a  range  of  frequencies  and  interpolating  back 
to  the  DC  value.  This  method  also  verified  the  value  for  the  DC  resistance.  For  each  of  these  scans  a  thin  plate 
of  aluminum  was  leveled  using  a  digital  protractor. 

A  series  of  air  scans  were  run  to  check  for  the  consistency  of  the  probe.  The  error  between  identical  air 
scans  taken  at  different  times  was  about  0.3%.  The  test  plate  was  then  placed  on  the  test  stand,  and  leveled 
using  the  same  digital  protractor  as  with  the  previous  test  plate.  The  probe  was  then  positioned  in  a  defect  free 
area,  with  the  center  of  the  probe  at  least  a  probe  radius  from  the  edge  of  the  notch.  The  probe  was  then  scanned 
parallel  to  the  notch,  using  the  same  scan  length  to  be  used  with  the  actual  notch.  The  purpose  of  this  was  to 
measure  subtle  variations  in  the  plate  conductivity  and  other  factors  that  may  impact  the  experiment.  Finally, 
the  scans  over  the  notch  were  performed.  These  were  run  such  that  the  probe  field  did  not  interact  with  the 
notch  at  each  end  of  the  scan. 

The  experiment  conducted  as  described  yielded  reproducible  measurements.  The  data  bracketing  process 
consisted  of  taking  two  separate  air  scans  before  running  a  test,  followed  by  anywhere  from  2  to  4  unflawed 
scans,  then  a  set  of  5  to  6  notch  scans;  this  was  then  repeated  in  reverse  order  to  verify  reproducibility. 

4.  Correction  for  Parasitic  Capacitance  and  Other  Non-ideal  Coil  Behavior 

A  final  step  needs  to  be  taken  after  data  collection.  Here  is  where  the  previous  measurements  taken  before 
the  scans  are  used.  This  parallel  network  comes  about  from  the  deviation  from  the  ideal  values  for  a  particular 
eddy  current  probe.  The  parallel  network  can  be  broken  down  into  three  main  components;  these  being  self¬ 
resistance  and  self-capacitance,  along  with  the  lead  capacitance.  This  correction  procedure  begins  with  using 
the  DC  resistance  and  inductance  values  to  find  ideal  admittance  [5]. 

y„=i/z„=i/(^,+j<»4)  (i) 

This  value  is  subtracted  from  the  air  admittance  value. 


Ya=UZa 


By  performing  this  subtraction,  the  admittance  of  the  parallel  network  is  obtained. 


(2) 


Yp  =  Ya-Y0  (3) 

The  admittance  of  the  probe  over  an  unflawed  region  of  the  plate  is  achieved  by  subtracting  the  admittance  of 
the  unflawed  region  from  the  parallel  network  admittance. 

Z£>„  =H(YU-YF)  (4) 

To  find  the  change  in  the  corrected  impedance,  this  value  for  the  corrected  unflawed  impedance  must  be 
subtracted  from  the  DC  value. 


A7  =7  -7 

CORR  ^ CORR  o 


(5) 


When  a  defect  is  seen  in  the  region  under  the  probe,  the  equation  below  gives  the  impedance  change. 


^DCORR=ll{YD-Yp)-ll{Yu-Yp) 


(6) 
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5.  Results 


In  this  section  the  results  from  the  experiments,  with  the  corrections  made,  as  well  as  the  simulations  from 
both  EC  SIM  and  VIC  3D  will  be  presented  for  two  of  the  plates.  In  looking  at  Figures  3a,  3c  and  3e  the  results 
have  a  good  agreement  in  shape  and  there  are  some  slight  discrepancies  in  the  magnitude.  The  majority  of  this 
discrepancy  is  a  shift  in  the  experimental  data,  either  up  or  down.  The  same  general  trends  can  be  seen  in 
Figures  3b,  3d,  and  3f.  Figure  3b  shows  a  legend  which  is  consistent  for  all  the  Figures  in  this  section. 


Figure  3a.  Resistance  for  6061  T6 


Figure  3b.  Resistance  for  7075  T6 


Figure  3c.  Reactance  for  6061  T6 


Figure  3d.  Reactance  for  7075  T6 


Figure  3e.  Impedance  for  6061  T6  Figure  3f.  Impedance  for  7075  T6 
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6.  Conclusion 


As  displayed  in  the  graphs  above,  there  is  qualitative  agreement  in  the  shape  of  each  of  the  curves.  There  is 
a  discrepancy  in  the  magnitude  of  the  plots  as  well  as  an  offset,  in  the  experimental  data.  This  can  be  seen  best 
in  Figure  3  a  and  Figure  3b.  One  would  expect  the  data  after  the  correction  to  approach  zero  away  from  the 
notch,  with  the  subtraction  of  the  unflawed  sections,  but  this  is  not  the  case.  This  could  indicate  that  the  probe 
interacted  with  the  edge  or  the  notch  at  ends  of  the  scan.  For  the  most  part  there  is  very  good  agreement  in  a  few 
of  the  plots,  in  both  magnitude  and  shape,  such  as  Figure  3c  and  Figure  3d  plate,  but  there  are  still  slight 
discrepancies  in  the  at  the  endpoints.  In  looking  at  Figures  3a,  and  more  importantly  3c  and  3e,  it  become  clear 
that  the  error  in  the  Figure  3e  is  due  to  the  shift  seen  in  Figure  3a. 

In  reviewing  the  results  of  the  two  simulation  codes,  some  interesting  trends  were  evident.  In  the  resistance 
component  shown  in  Figure  3a  and  Figure  3b,  there  is  good  correlation,  both  in  shape  and  magnitude.  This 
type  of  correlation  does  not  appear  in  the  other  plots.  In  Figures  3c-f  the  two  codes  produce  radically  different 
results.  In  comparing  the  simulation  codes  to  the  experimental  data  it  can  be  seen  that  in  the  majority  of  the 
plots,  the  results  generated  by  VIC-3D  are  more  in  agreement  with  the  experimental  data. 

One  outcome  of  this  investigation  was  an  appreciation  for  the  liftoff  control  necessary  for  this  type  of 
experiment.  The  original  study  made  use  of  a  spring-loaded  mechanism  to  keep  the  probe  at  desired  constant 
liftoff.  This  was  not  done  for  these  experiments.  In  addition,  the  mechanism  to  remove  any  tilt  form  the  probe, 
in  either  perpendicular  or  parallel  direction  to  the  scan,  was  rather  cumbersome  and  inaccurate,  which  made 
accurate  leveling  of  the  probe  somewhat  difficult. 
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